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ABSTRACT 


There are a number of studies in the literature where two 
dimensional finite difference analysis of complex groundwater 
basins are reported. However, not many studies on three 
dimensional models are available. In the present study a modular 
three dimensional finite difference model developed by the OS 
Geological Survey is implemented on a mini computer namely 
Micro VAX II system and also on DEC 1090, The model is then 
applied to a small basin in Anantapur district of Andhra Pradesh. 
The model being a sophisticated one, requires substantial amount 
of data for a few years. However, the data for only two years 
were available. Hence, the results are likely to be approximate. 
The main objective of the study was to demonstrate the 
applicability of three dimensional model to a groundwater basin. 

In addition to three dimensional finite difference model, the 
hydrological budget method is also used for the basin, as the 
basin is small. The advantage of the finite difference model over 
the hydrologic budget method is that it can take into account the 
spatial variation of hydrological and geological parameters and 
also complex boundary conditions. 

The results of the three dimensional finite difference model 
are voluminous. Only sample values of heads for layer number 
three at the end of December 1979 are reported in the thesis. 



CHAPTER I 


INTRODDCTION 


1.1 General: 

For the regions where the rainfall is very less, the 
availability of surface water becomes scarce. For such regions it 
is always advisable to study the possibility of groundwater 
resources and a proper development of the same should be 
considered. 

In planning and development of groundwater, care should be 
excercised to avoid undesirable effects of over-development. Some 
of the undesirable effects are: 

1. Gradual lowering of piezometric surface or groundwater table. 

2. Encroachment of saline water. 

3. Subsidence of land, etc. 

If groundwater is developed indiscriminately, it is possible 
that the amount of water withdrawn would be larger than the long 
term average groundwater recharged to the groundwater basin. This 
will result in a gradual decline of piezometric surface or 
groundwater table causing many shallow wells to go dry and making 
groundwater lift larger in deeper wells resulting in increased 
pumping cost. If the well fields are near sea coast, over 
development will cause encroachment of sea water into fresh water. 
Further, there is evidence of land subsidence in reaches where the 
aquifers are intersperced with compressible silty and clayey 
layers, from which large scale goundwater development has taken 
place. 
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Thus, whenever groundwat-er development, is considered, proper 
planning is required to avoid undesirable effects mentioned above. 
Such planning requires making suitable mathematical models of the 
basin and analysing the effects of development of groundwater. A 
large number of two dimensional mathematical models have been 
developed and applied to the groundwater basins. Relatively fewer 
number of three dimensional modeling methods have been reported in 
the literature. Recently, a general modular three dimensional 
finite difference groundwater flow model has been developed by the 
DS Geological Survey (1984). This modular program is quite 
versatile. In the present study this program was implemented on 
the Micro Vax II system. After fully implementing the program it 
was applied to one of the basins in the Anantapur district of 
Andhra Pradesh for which some data was available. 

1.2 Objectives of the sutdy: 

Given the data like rainfall, evaporation, canal flows, tank 
gauging observations, geological features and aquifer parameters, 
a model of groundwater basin can be constructed for evaluation of 
the groundwater potential and it can be simulated for various 
policy options and the spatial distribution of the heads can be 
predicted. 

The basin can be studied by simple hydrologic budgeting 
method for the basin as a single unit or by a mathematical model. 

4r 

The first method is very simple and hydrological balance was 
evaluated by considering simple analysis, where as in the detailed 
mathematical analysis a model of groundwater basin can be 
constructed regardless of its size and complexity. Even small 
details of geology and hydrology can be incorporated and studied 
by mathematical formulation. In constructing a groundwater basin 
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model, the basin is divided into smaller areas of rectangular or 
square blocks. The centre of the block is considered as a node for 
the mathematical treatment, assuming that the recharge or 
discharge into the model area occurs at that node. It is also 
assumed that the head in the block is same as that at the 
corresponding node. This gives the spatial variation (or 
distribution) of the groundwater - 

The main objectives of the study are: 

1 . Implementation of the DS Geological Survey modular three 
dimensional finite difference program. 

2. Application of hydrologic budget method for the basin. 

3. Application of the three dimensional modular program to 
the study area. 

1.3 Scope of the study: / 

The US Geological Survey model is a very general model with 
several options. It uses the programing language Fortran 66. It is 
applicable without much modifications on main frame computers such 
as CDC, DEC, CRAY etc, using Fortran 77 compilers. In the present 
study the above program was implemented on a mini computer namely 
Micro ¥ax II. Modifications necessary for this were incorporated 
in this program. 

Obviously for a three dimensional model the data requirement 
is much more comprehensive than for a two dimensional model. For 
example, we require not only storativity and transmissivity values, 
at different points in lateral space but also at differnet points 
in the vertical direction. Such detailed data are usually not 
available for any basin. In this sense the three dimensional 
model can only be applied to basins which are taken up for 
research study. It is however possible to use this program for two 
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dimensional problems also by using suitable options provided in 
the program. In the present study, therefore, the first objective 
was to implement the three dimensional program on the Micro Vax II 
and DEC 1090 the second objective was to apply this to a small 
catchment in Anantapur district of Andhra Pradesh. Further, the 
simple hydrological budget method is also applied for the basin. 

All the data that was available for this catchment was 
obtained. Based on the Lithology of the area, the data was 
collected from four different depths (layers). The length of the 
data was very short, only two years (1979,1980). Hence the results 
of application of the model will be approximate. The aim of the 
study is to demonstrate the applicability of the three dimensional 
model rather than to get reliable results for the basin which 
would not be possible because of the paucity of data. 

1.4 Organization of the study: 

This thesis has been divided into five chapters: 

% 

i. Chapter 1 gives a general intq^rduction about the objectives 

i 

and scope of the study. 

ii. Necessary review of literature is presented in Chapter 2. 

iii. Chapter 3 deals with the implementation of the program, 

iv. Application of a simple hydrologic model and the three 
dimensional finite difference groundwater flow model are 
discussed in Chapter 4. 

Results, conclusions and suggestions for further study are 
presented in Chapter 5. 


V. 



CHAPTER 2 


REVIEW OF LITERATURE 

A model is a tool designed to represent a simplified version 
of reality (Wang and Anderson, 1982). A mathematical/ model 
consists of a set of differential equations that are known to 
govern the flow of groundwater. The reliability of predictions 
from a groundwater model depends on how well the model 
approximates the field situation. 

Since no model can really reproduce all the variety in a real 
groundwater basin it is of some value to recognise how good the 
model may be for the purpose intended. The model itself can be 
used to reach a better understandng of the various aspects and 
phenomena involved in the question of model suitability 
(Thomas, 1973) . 

Several types of models can be used to study groundwater flow 
systems . They can be divided into three broad categories 
(Prickett, 1975): sand tank models, analog models (including 
viscous fluid models, and electrical models) and mathematical 
models (including analytical and numerical models). A sand tank 
model consists of a tank filled with an unconsolidated porous 
medium through which water is allowed to flow. A major drawback of 
sand tank model is the problem of scaling down a field situation 
to the dimension of a laboratory model. The observations measured 
at the scale of a tank model are often different from conditions 
observed in the field. Viscous fluid models known as the Hjile-Shaw 
or parallel plate models in which viscous fluid is made to flow 
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between two closely spaced parallel plates and the observations 
are made . Electrical analog model consists of boards wired with 
electrical networks of resistors and capacitors. They work 
according to the principle that the flow of groundwater is 
analogous to the flow of electricity. This analogy is expressed in 
mathematical similarity between Darcy’s law for groundwater flow 
and Ohm’s law for the flow of electricity. Changes in voltage in 
electrical models are analogous to changes in groundwater head. A 
drawback in electric analog model is that it is unique for each 

aquifer system. When a different aquifer is to be used, an entirly 

/ 

new electrical analog model must be built. A mathematical model 
consists of differential equations that are known to govern the 
flow of groundwater. 

Though sand tank models and analog models were extensively 
used prior to 1960 and ever since then, due to the wide 
availability of high speed digital computers, mathematical models 
have been favoured. The solution to the mathematical models can be 
obtained either by analytical or by numerical techniques. 
Unfortunately the equations of flow and continuity in the form of 
partial differential equations do not lend themselves easily to 
rigorous analytical solutions where boundaries are complicated. 
Further, the user must have a sufficient mathematical background 
to compute such solutions (Wang and Anderson, 1982). For example. 
Toth (1962) presented an analytical solution to the simple 
regional groundwater flow. The model consisted of a topographic 
high on one side and a valley bottom on the other, where water was 
discharging and the water table was gently sloping towards the 
valley. To obtain this analytical solution Toth (1962) introduced 
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several simplifying assumptions like the medium is an homogeneous, 
isotropic aquifer with linear water table configuration and he 
approximated the problem domain by a rectangle. None of these 
assumptions are necessary to obtain a nvimerical solution. To 
obtain solutions for such a simple problem the analytical method 
is complicated and for very complex groundwater problem it is not 
always possible to obtain a solution. To overcome this difficulty, 
one should resort to approximate methods using numerical 
techniques . 

At present there are two well known numerical methods which 
are in wide usage. They are the Finite Element Method (FEM) and 
the Finite Difference Method (FDM). One more method namely the 
Boundary Integral Equation Method (BIEM) is also being used. Of 
all the numerical methods, the finite difference method is the 
oldest and the most widely used technique. Ever since the 
classical paper by Richardson (1910), introducing finite 
difference as a method to calculate approximately the solution of 
partial differential equations was presented and later Shaw and 
Southwell (1941) applied this method for the solution of steady 
state seepage problem, this method has been extensively used in 
groundwater flow. 

The partial differential equation governing the flow of 
groundwater can be approximated numerically by finite difference 
technique or finite element method. In doing so, one replaces 
continuous variables with discrete variables that are defined at 
grid blocks (or nodes). Thus the continuous differentiaT equation 
defining hydraulic head everywhere in the aquifer is replaced by a 
finite number of algebraic equations that define hydraulic head at 
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specific points. The system of algebraic equations is generally 
solved using matrix techniques. This approach constitutes a 
numerical model and generally a computer program is written to 
solve the equations on a digital computer. 

Following Shah and Southwell (1941), Staltman (1956) and 
Fayers and Sheldon (1962) were among the early users of finite 
difference-based niimerical method for studying the regional 
groundwater basins. Since then, there have been many developments 
in numerical methods based on finite difference approach. Tyson 
and Weber (1964) used non-uniform polygonal discretisation of the 
aquifer area and the resulting finite-difference equations were 
solved by Gauss-Seidel method. In this study, a groundwater model 
was developed for a non-homogeneous but isotropic aquifer under 
confined or unconfined condition. The model was used for 
predicting the dynamic behaviour of a regional groundwater basin. 

Based on the changes of homogenity at different places of a 
basin, Fiering (1964) used a variable size rectangular grid while 
discretizing a nonhomogeneous and isotropic aquifer in both 
confined and unconfined situation. The set of finite difference 
equations was solved by an iterative technique. The model was 
aimed at predicting drawdown in the aquifer under variable well 
pumping rates and recharges. By a similar kind of model in which 
Gauss-Seidel iterative technique was used to solve the set of 
difference equations, Remson et.al. (1966) analysed the recharge 
effects of a proposed reservoir on a groundwater basin. 
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Peaceman and Rachford (1955) proposed a new method which has 
been found to result in significant advantages in solving 
groundwater flow problems. This was named as Alternate Direction 
Implicit (ADI) Method and in the later years the iterative ADI 
methods became very popular to solve the set of simultaneous 
finite-difference equations (Finder and Bredehoeft, 1968; Prickett 
and Lonnquist, 1968; Finder, 1970; Bredehoeft and Finder, 1970; 
Finder and Cooper, 1970; Rushton and Tomlinson, 1971; Bredehoeft 
and Finder, 1973; Lin, 1973). In these studies various types of 
groundwater problems such as transient flow in multiple layered 
aquifer involving storage in confining layers, effect of 
dispersion on steady movement of salt water front, and predicting 
future movements of contaminants were studied for various types of 
aquifer situations. 

A three dimensional model was described by Freeze (1971) 
which could handle heterogeneous, anisotropic, multilayer 
saturated-unsaturated and conf ined-unconf ined groundwater basins. 
The nodal sizes do not have to be constant but they must be block 
shaped (rectangular in two dimensions). This paper describes the 
use of the line successive over-relaxation technique to solve the 
problem numerically. Another interesting finite difference 
technique was discussed by Bibby and Sunada (1971) where the 
implicit central finite difference scheme was used and a solution 
was obtained by the Gaussian elimination procedure. This 
particular model was prepared to handle the leaky aquifer 


situation. 
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Any traditional approach may be appropriate for a system 
governed by a single partial differential equation but for a 
system comprising portions governed by effectively different 
equations may make the modeling process difficult. For any real 
groundwater system, it is likely that detailed modeling may 
require extensive computer capacity, CPU time and considerable 
amount of input data which consequently may prove to be an 
important restriction (Haddock, 1973). In particular, this 
difficulty stands for a large scale aquifer system where direct 
use of traditional technique may prove to be inadequate. The above 
difficulty is overcome by the method of decomposition (Dreizen and 
Haimes, 1977; Sharma and Lakshminarayana , 1986). In their study a 
large and complex aquifer system was decomposed into a number of 
subareas according to certain considerations involing 
geographical, geological and hydrological characteristics, 
administrative and operational judgements, or any other 
requirements associated with the particular need for the 
groundwater simulation model. 

Estimation of groundwater recharge is one of the important 
components of groundwater studies. Problems were encountered by 
Rushton et al. (1988), when trying to construct a recharge/ 
abstraction balance for the Liverpool sandstone aquifer. 
Potential recharge rates were calculated by a soil moisture 
balance method (Rushton and Ward, 1979), yet actual recharge to 
the aquifer had to be reduced by a factor which was dependent on 
the thickness and clay content of drift deposits lying between the 
soil and the aquifer. This was complicated further at grid points 
where the water table lay within the drift profile. This high- 
lights the need for comprehensive geological survey as well as 
water flux measurement to get a true picture of recharge to an 


aauif er . 
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The use of surface infiltration as a measure of recharge 
leads into with difficulties and uncertainities , partly because of 
the artificial effects of ponding, which is necessary for 
volumetric work and for ensuring potential rates of infiltration. 
Also, surface infiltration is generally more variable than 
transmissivity lower down the profile, which results from the 
nature of the surface, controlled as it is by the cumulative 
effect of raindrop and human compaction of the natural soil 
surface structure, and by the vegetation cover. 

A better approach involves estimation of soil water fluxes, 
allowing recharge to be estimated by the water balance equation. 
It is generally recognized that climate is one of the most 
important factors determining the amount of water loss by 
evapotranspiration from the crop. Apart from the climatic 
factors, evapotranspiration for a given crop is also determined by 
the type of crop and so are growth characteristics. Local 
environment, soil and soil water conditions, fertilizers, insect 
and disease infestations, agricultural and irrigation practices 
and other factors may also influence growth rates and result in 
evapotranspiration (FAO Irrigiation and Drainage Paper 24, 1975). 

Methods are used to predict evapotranspiration from climatic 
variables owing to the difficulty of obtaining accurate direct 
measurements under field conditions. Most prediction formulae use 
a differentiation between the components of climate and crop. 
Such formulae often have to be applied under climatic and 
agronomic conditions very different to those for which they were 
originally developed. 
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Evaporation pans provide a measurement of the integrated 
effect of radiation, wind, temperature and humidity on evaporation 
from an open water surface. In a similar fashion the plant 
responds to the same climatic variables but several major factors 
may produce significant differences in loss of water. 
Reflectivity of radiation from a water surface is only 5-8%, while 
that from most vegetative surfaces is 20-25% of solar radiation 
received. Daytime storage of heat within the pan can be 
appreciable and may cause almost equal distribution of evaporation 
between night and day, while most crops lose 95% or more of their 
24-hour loss during daytime hours. 

Reference crop evapotranspiration (ETo) can be predicted by 

ETo = Kp * Epan 

where ETo is in mm/day, Epan is pan evaporation in mm/day and 
represents the mean daily value of the period considered; Kp is 
the pan coefficient. 

The values of Kp are given in the FAO Irrigiation and 
Drainage Paper 24 (1975) for differnt conditions of humidity and 
wind, pan environment and type of pan and should be applied to 
pans located in an open environment with no crops taller than 1 m 
within some 50 m of the pan. Coefficients relating open water 
evaporation Eo to reference crop evapotranspiration ETo are 
presented in FAO Irrigiation and Drainage Paper 24 (1975). These 
values apply to shallow reservoirs and lakes with depths of less 
than 5 m and can be used to compute monthly Eo, once monthly ETo 
has been determined. The presented values apply equally to deep 
reservoirs and lakes in equitorial zones. 
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Lit.erat.ure review on some case studies are reported in the 
following section. 

Because of the importance of the integrated multi- 
disciplinary and basin wise groundwater potential, the CGWB {1983) 
has undertaken special projects based on socio-economic setup, the 
diverse hydrometeorological and hydrological settings for carrying 
out groundwater balance studies. The Noyil, Amaravati-Ponnani 
river basin projects of South India, for which extensive 
investigation was reported. The main source of recharge depends 
upon North-East monsoon where wide variation of monthly and annual 
rainfall and interannual variations are very high (ranging from 
20% to 90%). The hydrological investigation preceeded by well 
inventories (to a depth of 304m) indicated water bearing fractured 
zones. This type of report gives clear understanding of various 
aspects of groundwater balance evaluation and possible ways of 
improving the groundwater reserves. In this project the solution 
was arrived by using numerical techniques i.e. the 
finite-difference method. 

Beken and Byloos (1977) applied a monthly water balance model 
to Grote Nete test basin (553 sq km) in the North of Belgium. This 
low region has a complex geological structure, its boundaries are 
more or less unknown and deep infiltration into a deep aquifer is 
most likely to occur. This area is covered by several navigation 
canals. The inputs are monthly precipitation and Penman potential 
evapotranspi ration values. The model computes actual 
evapotranspiration, water storage in the basins, direct runoff and 
infiltration, base flow and total stream discharge, deep 
infiltration loss into the underlying aquifer and constant seepage 
from the canals. The author claims an accuracy of 0,2%. 



Schicht and Walton (1961) conducted studies for three basins 
in Illinois. Groundwater recharge runoff and evapotranspiration 
were determined for the study basins, from hydrologic and 
groundwater budgets using precipitation, stream flow and 
groundwater level data. The study basins are in North-Central, 
West-Southwestren and East-Central Illinois. Groundwater runoff 
under flood hydrogaphs, during rainless periods, groundwater 
evapotranspiration were determined from rating curves of mean 
groundwater stage versus groundwater runoff. Groundwater storage 
changes were computed from data on changes in mean groundwater 
stage by using a graph showing the relationship between gravity 
yields and average time of drainage. The gravity yield-average 
time of drainage graphs were prepared from the results of 
hydrologic budget studies made for periods during winter and early 
spring months when evapotranspiration and soil moisture changes 
were very small. Actual evapotranspiration is appraised from 
annual hydrologic budgets and were compared to mean annual 
potential evapotranspiration and annual water laws. 

The Central Ground Water Board (Southern Region, 1988) has 
drilled exploratory bore holes at different places in Ananthapur 
district of Andhra Pradesh and presented the lithology of the 
area. They have also conducted air tests and pumping tests to 
determine the aquifer parameters, also the drilling time logs of 
exploratory bore holes were reported as Basic Data Reports. These 
reports give some useful data necessary for the present study. 



CHAPTER 3 


IMPLEMENTATION OF MODULAR THREE DIMENSIONAL 
GI^UNDNATER FLOW MODEL OF THE D.S. GEOLOGICAL SURVEY 


3.1 MATHEMATICAL tK)DEL: 


The movement of groundwater of constant density through 
porous earth material may be described by the partial -differential 
equation 

dx J ^ &y [*^yy dy J ^ Wz~\^zz da J ~ ^s at (3-1) 

Where 


X, y and a are cartesian coordinates aligned along the major 

axes of hydraulic conductivity K , K , K ; 

x,x. yy z z 

h is the pieaometric head (L); 

W is a volumetric flux per unit volume and represents 
sources and /or sinks of water (T ^); 

S is the specific storage of the porous material (L ^); and 
t is time (T). 

In general, S„, K , may be functions of space 

s XX yy zz 

{S„=S„(x,y,z), and K ^=K„„(x,y,z), etc.) and h and W may be 
functions of space and time (h=h(x,y, z,t) , W=W(x,y,z,t) so that 
equation 3.1 describes groundwater flow under nonequilibrim 
conditions in a heterogenous and aniostropic medium. 
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In a finite difference model the domain is divided into a 
number of rectangular or square cells. Equation 3.1 is applied to 
each cell either by using finite differences for each differential 
term or by using a physical approach as below. 

Ah 

= ^s zrr ^3-2) 

where q^^is a flow rate into the cell (L^ t~^); 

Sg is the specific storage defined as the ratio of the volume of 
water which can be injected per unit volume of aquifer material 
per unit change in head {L~^) 

3 

AV is the volume of the cell {L ); and Ah is the change in 
head over a time interval, At. The term on the right hand side is 
equivalent to the volume of water taken into storage over a time 
interval At given a change in head of Ah. 

Figure 3.1 depicts a cell i,j,k and six adjacent aquifer 
cells i-l,j,k: i+l,j,k; i,j+l,k; i,j-l,k; l,j,k-l, and l,j,k+l. 
Flow into cell i,j,k in the row direction from cell 1, j-l,k (Fig. 
3.1), according to Darcy’s law, is given by 

^i. j-l,k‘^i, j,k 

"^i, j-l/2,k'" *®i, j-l/2,k^''i^^k 
where 

^i j-1/2 k volumetric fluid discharge through the face- 

' 3 -1 

between cells i, j,k and i, j'~l,k(Ii t )j j~l/2 k the 

hydraulic conductivity along the row between nodes i, j,k and 
i. j-l,k(L). 
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The index j-1/2 is used 1:0 indicate the space between nodes 
(Fig. 3.2). It does not indicate a point exactly halfway between 
nodes. For example, KR^ j-1/2 k ^®P^®sents hydraulic conductivity 
in the entire region between nodes i, j,k and i,j-l,k. 


Similar expressions can be written forapproximately the flow 
into or out of the cell through the remaining five faces, as: 


^i,d+l,k-^i, j,k 

%.j+l/2,k " ^i, j+l/2,k^°i^^k Arj^^/2 

^ i+l,j,k^ i,j,k 
'^i+l/2.j,k ^^i+1/2, j.k^^j^^k Ac^^^/2 ” 

^i-1, j,k-^i, j.k 

%-l/2,j.k ■ *^^i-l/2. j.k^^j^'^k 

^i, j,k+l-\. j,k 

*^i,j,k+l/2 = ^^i, j,k-H/2^^/'^i 

^i, j.k-l-^i, j,k 

%.j,k-l/2 " *^^i. j,k-l/2^^j^®i AVj^_i/2 


(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 


Each of the equations 3.3 -3.8 expresses flow through a face 
of cell i,j,k in terms of heads, grid dimensions, and hydraulic 
conductivity. Grid dimensions and hydraulic conductivity remain 
constant throughout the solution process so that the notation can 
be simplified by combining the constants into a single constant,' 
which multiplies head, called the "hydraulic conductance" or, more 
simply, the "conductance" , For example 

^^i , j -1/2 , k“ *®i , j-1/2 , k^^i^^k^^^ j -1/2 


(3.9) 



Where 
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j-l/2,k conductance in row i and layer k between nodes 
i,j-l,k and i,j,k (L^t~^). 

Conductance is the product of hydraulic conductivity and 
cross-sectional area of flow divided by the length of the flow 
path; in this case, the distance between the nodes. 


Substituting this expression into equations 3.3 to 3.8 yields 
j-l/2,k=^^i. j-l/2.k^^i, j-l.k'^i, j.k^ (3.10) 


‘^i.j+l/2,k= ^^i, j-i-l/2,k^^i. j+l,k ^i.j.k^ 


(3.11) 


‘^i-1/2, j,k= ‘^^i-1/2, j,k^^i-l, j,k' ^i, j,k^ 


(3.12) 


%+l/2,j,k" ^*^1+1/2, j,k^^i+l,j,k" ^i, j,k^ 


(3.13) 


‘^i,j,k-l/2" ^^i, j,k-l/2^^i, j,k-r ^i.j.k^ 


(3.14) 


%,j.k+l/2" ^^i, j,k+l/2^^i, j.k-t-l ^i. j.k^ 


(3.15) 


The above equations represent flow 
external flow into or out of the cells, 
well pumpage can be expressed as 

®i,j,k,n~ ^i, J,k,n\, j,k'^ %,j,k,n 


between the cells. The 
such as river recharge or 

(3.16) 


Where 

a. i V represents flow from the n-th external source into 

X I j j ic I n 

cell i, j,k(L^t"^), and p^ j,k,n %, j,k,n constants 

(L^t ^ and L^t respectively) . 
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Substitution of equtions 3.10 to 3.16 into equation 3.2 yields, 

‘^‘^i-l/2,j.k‘*'i-l.j.k'*‘i,j,k> ■"°°l+l/2,J,k‘’'l+l,J,k‘*'i, j,k> ■" 

‘^''i,j,k-l/2<‘'i,j.k-l'l'i.j,k> ■"‘^''i,j,k+l/2<*‘i,J,k+l‘J'i, j,k^ ^ 




(3.17) 


where jj. is the sum of all j jj. from different external 
sources and j jj. is the sum of all components j ^ Using 
backward differencing for time, we get 


^®i , j-1/2, k^^i , j-1 , k"^i , j , k^‘^‘^®i , j+1/2, k^^i , j+1 . k ^i . j . k)"^ 

^^i-1/2, j,k^^i-l, i,k~^i, j,k^ ''’^^i+1/2, J,k^^i+1, j,k"^i, j,k^ 
^^i, j,k-l/2^^i, j,k-l~^i, j,k^ j,k+l/2^^i, j,k+l"^i, j,k^ 


P. . » hj j ^4 < V 

XitJjrJS X|J||>K jL p Ji $ *k 




.,m % 

^^i, j,k'^i, j,k^ 


(3.18) 


t -t ^ 
m m-1 


where h 


i, j,k 


is the head at the end of time interval m and 


h i j k known head at the end of the previous time 

interval m-1. For example, if we apply equation 3.18 for time 
interval 2, we get 

^^i, j-l/2.k^^i, j-l,k'^L j,k^'^^®i, J+l/2,k^^i, j+l,k'^L 

‘^^i-1/2, j,k^^Ll, j.k'^L j,k^ ■^^^i+1/2. j,k^^i+l, j,k"^i, J,k^ ^ 

^’^i, j,k-l/2^^i, J,k“l"^i. j,k+l/2^^i, J,k+l’'^i, j,k^ 


^i. J .k\. j.k"^ ^i.j,k 




^^i,J.k ^i, j,k^ 


■^2"’*'l 


(3.19) 
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Equation 3.18 can be written in the following form 


CV 


i, j,k-l/2^i, j,k-l ^^^i-1/2, j,k^i-l, j,k ''’^^i, j-l/2,k^i, j-l,k 


-CC. 


-CR, 


j.k-1/2 ""i-l/2,j,k ^^i , j-l/2,k‘‘^^i, j+l/2,k 

CC^ . , ^ ,.-CV^ . ,„+ HCOF^ . ,_)h“ j ,_+CR^ ^ /r, ,_h“ 


'i+l/2.j,k '"’i.j,k+l/2 


+CC 


.m 


i,l,k'“i, j,k ""“i, j+l/2,k“i, j+l,k 


+CV 




i+1/2. j.k“i+l. j.k^^''i. j.k+l/2"i, j,k+l = ^^i.j.k 


(3.20) 


where 

HOOF 


RHS 


i, jk 
i, j,k 


^i,j,k ^^^i, 


(lV^) 


and (lVI) 




(L^) 


(3.21) 


Equation 3.20 is applied to each cell in the domain. We then get a 
set of linear algebraic simultaneous equations which should be 
solved at each time step. These simultaneous equations are usually 
solved by iterative methods. 

3.2 FBOGRAH DESIGN: 

The computer progreun to simulate Three Dimensional Finite 
Difference Groundwater model was developed by the Water Resources 
Divison, Ministry of Interior of the United States Geological 
Survey (Me Donald and Harabaugh, 1984). The following description 
closely follows the above reference. This model combines all the 
requirements for the evaluation of groundwater flow. The 
components like rainfall recharge, river recharge/discharge, well 
recharge /discharge, groundwater evapotranspiration, general head 
flow boundaries are all included in the model. Provison has also 
been made for inclusion of several other new options without 
affecting the overall structure of the program. The basic structure 
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of the program consists of a main prograim and a number of highly 
independent modules which perform the individual tasks for which 
they are meant. The main program invokes a particular module based 
on its necessity to compute the requirements provided in that 
module. The unwanted modules will not be utilized while simulating 
the model. Figure 3.3 shows the organization of the simulation. 

The overall program structure with brief explanation is given 
in figure 3.4. The main program consists of a number of procedures 
namely DEFINE, ALLOCATE, READ and PREPARE, etc. Each simulation 
has three nested loops : a STRESS PERIOD LOOP, within which there 
is a TIME STEP LOOP and which inturn consists of an ITERATION 
LOOP. 


In the DEFINE Procedure, the problem to be simulated is 
defined specifying the size of the model, the type of simulation 
(transient or steady state), the number of stress periods, 
hydrologic options and the solution scheme desired. In the 
ALLOCATE procedure, memory space required by the program is 
allocated. In the READ AND PREPARE procedure, all the data that 
are not functions of time are read including boundary conditions, 
initial heads, transmissivity /hydraulic conductivity, specific 
yield/storage coefficient, elevation of layer top and bottom and 
parameters of the specified solution scheme. Each package 
consists of all modules associated with a particular hydrologic 
feature and a solution method. For example, the RIVER package, 
performs all the tasks needed for simulation of the river flow 
effect on the model, similarly other packages like WELL package, 
DRAIN package, RECHARGE package, etc., perform their respective 


tasks . 
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3.3 IMPIJMENTATI(»I OF THE PBOGStAM : 

The computer program originally developed by US Geological 
Survey, Water Resources divison was written in Fortran 66. It is 
implemented in the Micro Vax II system which has a fortran 77 
compiler and also DEC - 1090 at IIT Kanpur. Some of the 
modifications made while implementing the program are given below. 

1) The program listing given does not contain any file "OPEN" 
statement. The Micro VAX II system resulted • incertain problems 
during the execution and debugging of the program and hence the 
same is included. 

OPEN(UNIT= 1,FILE=’ INBAS.DAT’ , STATUS= ’ OLD ’ ) 
0PEN{UNIT=16,FILE=’ IOUT.DAT’ , STATUS=’NEW’ ) 

2) As the accuracy of the model till sixteen digits is not 
necessary, the DOUBLE PRECISION statements are deleted so as to 
save memory and also to reduce the computational time 
considerably. To effect this change the stringe is modified as, 

NLAY * NROW * NCOL 

3) Data statements in Micro VAX II system do not 
require ( ’ ) before and after end of data string and the same was 

modified. 

The program as modified above was compiled without any 
compilation errors. 

3.4 VALIDATION OF THE PROGRAM: 

The users manual contains one set of test data with typical 
output. The program was executed with the same sample input data 
and the results obtained are the same as given in the manual. 

A summary of the results is presented in Table 3.1 for comparison. 
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Table 3,1 Comparison of Eesults of Test- Data 

lOL I BIETRIC £IIDfflET IQE ENTIRE MODEL AI EM QE TIME STEP 1 

m STRESS EERIOD 1 


REPORTED Micro VAX II DEC 1090 

IN THE MANUAL 


CUMULATIVE VOLUMES L ** 3 


mi 


STORAGE 

= 0.00000E+00 

0 , 00000E+00 

0.00000E+00 

CONSTANT HEAD 

= 0.00000E+00 

0.00000E+00 

0.00000E+00 

WELLS 

= 0.00000E+00 

0.00000E+00 

0.00000E+00 

DRAINS 

= 0.00000E+00 

0.00000E+00 

0.00000E+00 

RECHARGE 

= 0.13608E+08 

0.13608E+08 

0. 13608E+08 

TOTAL IN 

= 0.13608E-t-08 

0.13608E+08 

0.13608E+08 

QU-Il 




STORAGE 

= 0..00000E+00 

0.00000E+00 

0.00000E+00 

CONSTANT HEAD 

= 0.43265E+07 

0.43265E+07 

0.43265E+07 

WELLS 

= 0.64800E+07 

0.64800E+07 

0.64800E+07 

DRAINS 

= 0.28011E+07 

0.28011E+07 

0.28011E+07 

RECHARGE 

= 0.00000E+00 

0.00000E+00 

0.00000E+00 

TOTAL OUT 

= 0.13608E+08 

0.13608E+08 

0.13608E+08 

IN - OUT 

= 184.00 

395.00 

403.63 

PERCENT DISCREPANCY = 0.00 

0.00 

0.00 

RATES EQE THIS UM SIEE L ** 3/T 



m 




STORAGE 

= 0.00000E+00 

0.00000E+00 

0.00000E+00 

CONSTANT HEAD 

= 0.00000E+00 

0.00000E+00 

0.00000E+00 

WELLS 

= 0.00000E+00 

0 . 00000E+00 

0.00000E+00 

DRAINS 

= 0.00000E+00 

0.00000E+00 

0.00000E+00 

RECHARGE 

= 157.50 

157.50 

157.50 

TOTAL IN 

= 157.50 

167.50 

157.50 

OUT- 




STORAGE 

= 0.00000E+00 

0.00000E+00 

0.00000E+00 

CONSTANT HEAD 

= 50.075 

50.075 

50.075 

WELLS 

= 75.000 

75.000 

75.000 ■ 

DRAINS 

= 32.420 

32.420 

32.420 

RECHARGE 

= 0.00000E+00 

0.00000E+00 

0.00000E+00 

TOTAL OUT 

= 157.50 

157.50 

157.50 

IN - OUT 

= 0.21210E-02 

0.045776E+-02 

0.46711E-02 


0.00 0.00 0.00 


PERCENT DISCREPANCY 
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Figure 33 -Cell i.j.k and indices for the six adjacent cells. 
(Adopted irom He Donald and Harabaugh, 1984) 
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Figure 3.3 /-Primary modules organized by procedure and 

package. 
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DEFINE — Read data specifying numt)ef of rows, 
columns, layers, stress periods, and maiof program 
options 

ALLOCATE — Allocate space m the compute! to 
store data. 

READ AND PREPARE — * Read data which is constant 
throughout the simulation Prepare the data by 
performing whatever calculations can be made at 
this stage 

STRESS — Determine the length of a stress period 
and calculate terms to divide stress periods into time 
steps. 

READ AND PREPARE — * Read data which changes 
from one stress period to the next. Prepare the data 
by performing whatever calculations can be made at 
this stage. 

ADVANCE — Calculate length of time step and sat 
beads at beginning of a new timt step actual to heads 
calculated lor the end of the previous time step. 

formulate — Calculate the coefficients of the 
fmiie difference equations for each cell. 

APROXIMATE — Make one cut at approximating a 
solution to the system of (intte dtfierence equations. 

OUTPUT CONTROL Determine aether results 
should be written or saved on disk for this time step. 
Send signals to the BUDGET and OUTPUT pro- 
cedures to indicate exactly what information should 
foe put out 

BUDGET — Calculate terms for the overafi vofu- 
metric budget and calculate and save cell-foy-cell 
flow terms for etch component of Bow. 

OUTPUT — Print and save heads, drawdown and 
overall volumetric budgets in accordance with 
signals from OUTPUT CONTROL procedure. 


Figure 3.4 -Overall program structure. 

(Adopted from Me Donald and Harabaugh, 1984) 













CHAPTER 4 

APPLIGATKMI OF THE MODEL 

4.1 Study Area 

Th^ Hiociel described in the previous chapter is applied to a 
small catchment (or basin) namely Havaligi catchment in Anantapur 
District of Andhra Pradesh. The Havaligi catchment lies between 
latitude 14^ 50’ - 15° 05’ and longitude 77° 05’ - 77° 20’ 

5 > 

covering an area of 267 km . 

The data for the above catchment were collected from Andhra 
Pradesh State Irrigation Development Corporation and Andhra 
Pradesh Ground Water Department. 

Generally, for any groundwater studies we need data for 
atleast about five to ten years. However, in the present case 
data were available for only two years (1979 - 1980). The types 
of data collected are: 

1 . Topography of the area 

2. Daily Rainfall 

3. Daily evaporation (Pan) 

4. Daily tank water levels and outflows 

5. Canal length, width and its discharge 

6. Borehole particulars 

a) Lithology 

b) Static water levels 

c) Diameter of bore and depth drilled 

7. Storage coefficient and transmissivity values 

8. Command area particulars 



Though we do not have agricultural data for the period of 
study, suitable assumptions based on present agricultural practices 
were made. 

4.1.1 Topography: The catchment areas of Havaligi and Urvakonda 
tanks (anicuts) are marked on the toposheets (Fig 4.1), 
delineating the catchment by joining the ridge lines. The streams 
and tanks present within each catchment area are also marked. The 
general topography can be easily seen from the Fig. 4,1. Basically 
the catchment under study is plain with very flat slope. 

4.1.2 Rainfall Data: There are fifteen raingauge stations in the 
Havaligi catchment. The rainfall data for the stations is given 
in Table 4.1. Though daily rainfall data for two years is avai- 
lable, only monthly values are given in Table 4.1. Some missing 
rainfall data were estimated by using the weighted average method. 

4.1.3 Pan Evaporation Data: The daily pan evaporation data were 
recorded at the hydrological station in the basin using ISI pan 
evaporimeter which is similar to D.S. Class A pan both in 
construction and in dimensions. These values are represented by a 
graph ( Fig .4.2). 

4.1.4 Tank Water Levels: There are two tanks in the catchment 
area. Daily water levels of the tanks are available and the 
contours of the tanks are also available, from which we can 
calculate the water spread area, which is useful for calculating 
evaporation losses from the tanks and seepage losses. Also, we 
have the daily outflows from the tanks through the sluices and 
over the surplus weirs. Figure 4.3 shows the contour map of 



30 


Havaligi tank and figure 4.4 slw- .iap of the 
Urvakonda tank. Figures 4.5 and 4.6 respectively show the stage 
versus water spread area and stage versus cumulative volumes of 
Havaligi and Urvakonda tanks. 

4.1.5 Canal Particulars: Two major unlined canals pass through 
Havaligi catchment. These are Tungabhadra High Level Canal and 
Guntakal Branch Canal. These are shown in Fig. 4.1. For water 
balance study one would require the amount of water that enters 
the basin through canal and the amount of water that leaves the 
basin through the same canal , This would enable computation of 
canal losses over the basin. However, such information is not 
available for the present basin. The seepage losses from the 
canals are calculated by using Darcy’s law and by making 
assumptions about canal bed material. 

4.1.6 Borehole Data: Data on several boreholes in the catchment 
area were available. The lithological data collected while 
drilling wells are available. A typical lithology is given in 
Fig. 4.7. 

4.1.7 Command Area Particulars: The command area of Urvakonda 
tank is presented in fig 4.8. The type of crops grown in each 
season were not available. However, assumptions have been made 
about these based on the crops grown in the basin presently. 
Table 4.2 shows the type of crops and their crop coefficients. 

A groundwater balance study is required for evaluating' 
groundwater resources. In the following sections two methods are 
presented. In the first method, the overall hydrological balance 
for the basin as a whole including both surface water and 
groundwater is presented. In this method, the spatial variation 
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within the basin in hydrological parameters as well as aquifer 
parameters are not taken into account. In the second method, a 
detailed three dimensional model is used which takes into account 
all the variations in the basin ( Lakshminarayana 1989). " 

4.2 Hydrologic Budget Model 

Groundwater is a resource which is being developed more and 
more. Groundwater storage, in contrast with surface storage, is 
unseen. There are several questions one has to consider in 
connection with groundwater development: How much water can be 
drawn safely? What happens if more water than this is drawn? How 
can we increase the supply to groundwater basin? 

4.2.1 Hydrologic Balance 

In order to determine how much water can be taken out of 
groundwater safely without causing undesirable effects, one has to 
do groundwater balance or groundwater budget study. A groundwater 
basin is a part of the total hydrologic basin which includes both 
surface water and subsurface water. The general form of the 
hydrologic balance includes a statement of all the waters entering 
and leaving the basin. It may be expressed as 

Inflow - Outflow = Change in storage 

The inflow includes all components such as precipitation, 
surface inflow, subsurface inflow, water imported from outside 
the basin, etc. The outflow includes all items such as surface 
outflow, subsurface outflow, consumptive use, water exported 
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outside the basin, etc. The change in storage includes change in 
surface water storage and change in subsurface storage. The change 
in subsurface storage should account for both soil moisture 
storage and groundwater storage. 

Because of interrelation of surface and groundwater, all 
inflows and outflows of the basin must be included in the balance 
equation. Often certain items can be omitted because they are 
negligible. Any consistent unit of volume and time can be 
employed, say million cubic meters per year. The surface inflow, 
outflow, precipitation, imported and exported water are all 
measurable. The consumptive use or evapotranspiration can be 
estimated by determining the various types of crops grown and 
other types of vegetation such as forests in the basin. The sum of 
the products for each land use, area times consumptive use yields 
the total consumptive use. 

The changes in the surface storage can be evaluated directly. 
The subsurface inflow and outflow are diffcult to estimate. If it 
is possible to choose the basin such that no subsurface flow 
crosses into and out of the groundwater basin then these terms are 
zero. If sufficient data are available it would be possible to 
estimate them by Darcy’s law. 

By selecting time intervals in which the amount of water in 

) 

unsaturated zone is nearly equal, only changes in saturated zone 
need be considered. From measurement of water level in wells, maps 
of groundwater changes can be prepared. The product of change in 
water level times specific yield times area gives change in 


groundwater storage. 
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4.2.2 Groundwater Balance 

A specialized form of hydrologic balance equation in which 
all items pertaining to a groundwater reservoir is called a 
groundwater balance or budget equation. Recharge is the main 
inflow component for a groundwater basin. If there are inflows and 
outflows across the boundary they can be estimated by using 
Darcy’s law. The main outflows are groundwater runoff, groundwater 
evapotranspiration and groundwater pumpage, A groundwater balance 
equation can be written as 

P - R -ET - W= AH. S 

0 g y 

where P = groundwater recharge 
R = groundwater runoff 
W = groundwater pumpage (withdrawal) 

ET = groundwater evapotranspiration 

W 

AH = change in groundwater level 
Sy = Specific yield 

Again consistent system of units should be used in the above 
equation. 

4.2.3 Groundwater Recharge 

The main source of groundwater is precipitation which may 
infiltrate into the ground and finally reach groundwater table, or 
may enter surface streams and seep from these channels to the 
groundwater. Interception by vegetation and soil moisture must be 
satisfied before any large amount of rainwater can percolate into 
groundwater. Only prolonged period of heavy precipitation can 
supply large quantites of water for groundwater recharge . 



Groundwater recharge is an intermittent and irregular process. 
Water used in irrigation in excess of requirement of crops may 
also become groundwater recharge. Groundwater levels are 
usually at their lowest in the month of May i.e., just before 
monsoon and will reach their highest level in the month of 
October, just after the monsoon in this country. Groundwater 
development must be made in a planned way. On a long term basis 
the amount of groundwater taken from a groundwater basin should 
not be more than the natural groundwater recharge reaching the 
aquifer. 

4.2.4 Groundwater Evapotranspiration 

In areas where water table is shallow and the capillary 
fringe extends to the land surface, large quantites of groundwater 
are discharged into atmosphere by evaporation and transpiration by 
plants. This process is called evapotranspiration. 
Evapotranspiration depends primarily on meteorological factors, 
available soil moisture and groundwater , and type of soil and 
vegetation. At many places, evapotranspiration is small during 
winter months and reaches a maximum during summer months. In 
summer months evapotranspiration is very effective in reducing 
groundwater runoff. With the same mean groundwater stage 
groundwater runoff is much less in summer than in winter. Separate 
rating curves, that is groundwater stage versus groundwater runoff 
curves, must be prepared for summer and winter months. The 
difference in groundwater runoff between the two curves can often 
be used as a rough estimate of groundwater evapotranspiration. 

The following inputs and outputs are computed for applying 
the hydrologic budget equation to the Havaligi catchment. 
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4.2.5 Inputs 


4.2.5a Rainfall: Since there is variation in the distribution of 
rainfall over the basin as seen from the rainfall data recorded at 
different raingauge stations, the basin is divided into a a number 
of Thiessen polygons (Fig 4.9) and the rainfall within each of the 
Thiessen polygons is assumed to be the same as recorded by the 
raingauge located within the polygon. This method is justified 
because the basin is not located in the mountaneous region and the 
areas of each of the Thiessen polygons is small ( 3 % to 14% of the 
total catchment area). The total volume of the rainfall 
contribution is determined by 


I 


P 


^ 1^1 * ^ 2^2 ^ 


+ A P 


n n 


These results are shown in Table. 4.3 and the area of each of 
the Thiessen polygon are shown in Table 4.4 


4.2.5b Seepage from Canals and Tanks: The wetted area of the 
tanks are taken as the areas of water spread. For the tanks, the 
water spread area is computed as the mean of maximum and minimum 
water spread areas in each month. The water spread areas are 
obtained from the stage vs. water spread area curves (Fig 4.5 and 
4.6). Assuming the vertical hydraulic conductivity of the tank bed 
and the thickness of the bed material, the recharge to the 
groundwater is evaluated, using Darcy’s law. They are shown in the 


Table 4.5. 
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4.2.5. C Surface and subsurface inflows: If "there are any surface 
and subsurface inflows in"bo "the s"budy area, "they are also "to be 
considered for evaluation. The only surface inflow in"to "the s"budy 
area is "through canal and i"t is already considered in "the previous 
section. Due to lack of sufficient information about the 
subsurface inflows it is assumed that the subsurface boundary is 
congruent with the surface boundary of the area. Hence the 
subsurface inflows are absent. 

4.2.6. Outputs 

Not all the precipitated water reaches the groundwater. There 
are losses like Evaporation, Evapotranspirat ion, surface runoff 
etc. 


4.2.6a Evaporation and Evapotranspiration: Based on the pan 
Evaporation values, the reference evapotranspiration and potential 
evapotranspiration can be calculated 

ET = E * K 
0 ^ 

PET = ET^ * K 
0 c 

where K is pan coefficient. 

P 

K is crop coefficient, 
c 

ET 0 is reference evapotranspiration. 

PET is potential evapotranspiration. 

E is pan evaporation. ■ 

The values of K and K are taken from FAO Irrigation and 

P C — , ■ 

Drainage paper 24 (1975) and FAO Irrigation and Drainage paper 33 
(1978) and the calculated values of the PET are given in 


Table 4.6. 
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4.2.6b Surface Runoff: The surface outflows from the catchment 
area are available and are shown in Table 4.7. 

4.2.7 The Hydrologic Budget: 

Based upon the above mentioned inputs and outputs for each of 
the month, the yearly hydrologic budget is made to determine the 
groundwater recharge. Due to lack of information about the soil 
moisture, the same was not included in the balance equation. 
However, this will not affect the yearly hydrologic budget, since 
it is assumed that the soil moisture will be same for the same 
seasons of different years. The results are discussed in Chapter 5. 

4 . 3 The Three DiBtensional Finite Difference Model : 

In the previous section we have used hydrologic balance for 
the basin as a whole and obtained the toal volumetric budget per 
annum. In the present section we will include the variations 
within the basin form point to point. For this it is necessary to 
use three dimensional finite differnece approach. The basic 
equations for this method have already been given in Chapter 3. 

4.3.1 Discretization of the basin: 

In order to apply the finite-difference method, the basin was 
discretized into rows, columns and layers. The basin under present 
study is Havaligi basin. This basin has been discrftized into 31 
rows, 40 coltnnns and 4 layers (Fig 4.10). The space between the 
rows and columns is same forming a square block of each layer 
which is same throughout that layer and it differs form other 
layers. These layers are the actual thicknesses of each of the 
geological formation present in the basin, whch was taken from 


borehole data. 
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4.3.2 Boundary Condi'bions: 

The boundary conditions are applied to the basin. The 
boundary conditions considered are of three types and are 
specified by codes which were referred as 

IBOUND < 0 constant head cell. 

IBOUND = 0 inactive cell. 

IBOUND > 0 variable head cell. 

The constant head cell is one in which the variation of head 
is absent. And the variable head cell is the one in which the 
groundwater head changes in conduction with the various inputs, 
outputs and interaction with boundaries. The discretized basin 
with boundary conditions are shown in fig 4.11. 

4.3.3 Initial Ck>nditions 

Though data for each cell is not available, but sufficient 
number of well data could over come this difficulty. Based on the 
available information of static water levels in the wells, the 
water table contours were drawn fig 4.12 and the same data is 
given as input for initial heads in each of the cells. 

4.3.4 Rainfall Recharge 

The net effective rainfall recharge is computed based on the 
availabe data like rainfall, evapotraspi ration and the surface 
runoff and the same is given as effective recharge to each of the 
cells. They are presented in Table 4.8. 
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4.3.5 CaBal and Tank. Beckarge : 

There are two canals passing through the basin under study 
and the cells through which they pass are shown in the Fig 4.1. The 
data required for the evaluation of canal seepage is shown in 
Table 4.9. Though the program does not have provision for 
computing the tank reacharges/discharges , we have utilized the 
canal package for this evaluation because the mathematical 
formulation for both cases is the same. 

4.3.6 Well Pumping: 

The location of wells and their particulars such as depth 
drilled, quantities of water pumped for each time period are 
given as the input for evaluation Table 4.10. 

4.4 Application of the model: 

For the detailed evaluation of the Havaligi basin using three 
dimensional finite difference method, the study area is 
discretized as explained earlier. Considering a time step of one 
month in each stress period the simulation is carried out using 
Slice Successive Over Relaxation (SSOR) solver package using the 
above mentioned data. The volumetric budget for each stress 
period and the head in each cell for each layer are estimated. The 
drawdown in each cell and in each layer are also obtained. The 
results are discussed in Chapter 5. 
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Table 4.2 Types of crops and their crop coefficients 


CROP INITIAL 

STAGE 

DEVELOPMENT 

STAGE 

MID 

SEASON 

LATE 

SEASON 

TOTAL 

PERIOD 

RICE 

1.15 

1.50 

1.30 

1.05 

1.20 

— >■ 

30 

30 

46 

55 

160 

GROUNDNUTS 

0.50 

0.80 

1.10 

0. 85 

0.80 

^ 

26 

35 

45 

25 

130 

MILLETS 

0.40 

0.75 

1.15 

0.80 

0.85 


16 

25 

40 

25 

105 

COTTON 

0.60 

0.80 

1.26 

0.90 

0.90 


30 

50 

60 

55 

196 


— ^ represent length of each growth period in days - 



TABLE 4.3 MONTHLY RAINFALL VOLUMES IN EACH THIESSAN POLYGON ( Million Cubic metrea 
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TABLE 4.6 


POTENTIAL 

EVAPOTRANSPIRATICMI 

VALDES 

Month 

1979 

1980 


(million cubic 

meters) 

January 

1.3073 

1.4488 

Febraury 

1.3302 

0.8983 

March 

0.0000 

0.0529 

April 

0.0000 

0.0000 

May 

0.0000 

0.0000 

June 

5.4554 

5.4505 

July 

6.9338 

5 . 6650 

August 

5.8690 

3 . 5010 

September 

7.3740 

9.6124 

October 

4.7453 

4.8105 

November 

1.8185 

0.8013 

December 

1.4499 

0.7270 



TABLE 4.7 


OUTFLOWS FROM HAVALIGI TAMK 


Month 

1979 

(million 

1980 

cubic meters) 

January 

0.0419 

0.0414 

Febraury 

0.0011 

0.0036 

March 

0.0000 

0.0000 

April 

0.0000 

0.0000 

May 

0.0000 

0.0000 

June 

0.0000 

0.0000 

July 

,0 . 0000 

0.0000 

August 

0.0692 

0.0736 

September 

0.0215 

0.0181 

October 

0.1428 

0.0000 

November 

0.1464 

0.0000 

December 

0.1219 

0.0000 



TABLE 4.8 EFFECTIVE RAINFALL RECHARGE TO THE GROUNDWATER 
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TABLE 4S CANAL PARTICULARS 

(JANUARY 1979) 


Layer 

Row 

Column 

Stage 

(m) 

Conductance 
( n^month) 

Bottom , 

Elevation (bgl) 

2 

19 

7 

1.7 

365.60 

3.0 

2 

19 

8 

1.7 

365.60 

3.0 

2 

18 

.9 

1.7 

365.60 

3.0 

2 

18 

10 

1.7 

365.60 

3.0 

2 

18 

11 

1.7 

365.60 

3.0 

2 

18 

12 

1.7 

365.60 

3.0 

2 

18 

13 

1.7 

365.60 

3.0 

2 

19 

14 

1.7 

365.60 

3.0 

2 

19 

15 

1.7 

365.60 

3.0 

2 

20 

16 

1.7 

365.60 

3.0 

2 

21 

17 

1.7 

365.60 

3.0 

2 

21 

18 

1.7 

365.60 

3.0 

2 

22 

19 

1.7 

365.60 

3.0 

2 

23 

20 

1.7 

365.60 

3.0 

2 

24 

21 

1.7 

365.60 

3.0 

2 

24 

22 

1.7 

365.60 

3.0 

2 

24 

23 

1.7 

365.60 

3.0 

2 

24 

24 

1.7 

365.60 

3.0 

2 

24 

25 

1.7 

365.60 

3.0 

2 

24 

26 

1.7 

365.60 

3.0 

2 

24 

27 

1.7 

365.60 

3.0 

2 

24 

28 

1.7 

365.60 

3.0 

2 

23 

21 

1.7 

268.70 

3.0 

2 

22 

22 

1.8 

268.70 

3.0 

2 

21 

22 

1.8 

268.70 

3.0 

2 

20 

22 

1.8 

268.70 

3.0 

2 

19 

22 

1.8 

268.70 

3.0 

2 

18 

22 

1.8 

268.70 

3.0 

2 

17 

23 

1.8 

268.70 

3.0 

2 

16 

23 

1.8 

268.70 

3.0 

2 

15 

23 

1.8 

268.70 

3.0 

2 

14 

22 

0.9 

234.20 

3.0 

2 

13 

22 

0.9 

234.20 

3.0 

2 

13 

21 

0.9 

234.20 

3.0 

2 

12 

20 

0.9 

234.20 

3 . 0 

2 

11 

19 

0.9 

234.20 

3.0 

2 

11 

18 

0.9 

234.20 

3.0 

2 

10 

18 

0.9 

234.20 

3.0 

2 

9 

17 

0.9 

234.20 

3.0 

2 

8 

17 

0.9 

234.20 

3.0 

2 

7 

17 

0.9 

234.20 

3.0 

2 

6 

17 

0.9 

234.20 

3.0 


^ bgl below ground level 
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Table 4.10 Wells Location and Discharges 
for January 1979 


JANUARY 1979 


LAYER 

ROW 

COLUMN 

DISCHARGE (m/day) 

3 

2 

11 

50.0 

3 

5 

12 

37.5 

3 

6 

27 

50.0 

3 

7 

8 

60.0 

3 

8 

15 

45.0 

3 

8 

20 

60.0 

3 

8 

33 

50.0 

3 

9 

24 

50.0 

3 

9 

29 

50.0 

3 

10 

38 

50.0 

3 

11 

17 

52.5 

3 

11 

26 

50.0 

3 

12 

6 

50.0 

3 

14 

23 

50.0 

3 

14 

37 

50.0 

3 

15 

12 

75.0 

3 

16 

15 

60.0 

3 

16 

20 

50.0 

3 

16 

24 

50.0 

3 

16 

32 

50.0 

3 

17 

35 

50.0 

3 

18 

19 

50.0 

3 

18 

29 

75.0 

3 

20 

13 

50.0 

3 

20 

24 

90.0 

3 

21 

28 

50.0 

3 

22 

15 

50.0 

3 

23 

24 

90.0 

3 

24 

19 

50.0 

3 

25 

22 

50.0 

3 

26 

26 

50.0 
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Pl^. 4.1 Topography of Havallgi Catchment Area 
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(Stage Vs area and Stage Vs Cumulative Volume) 
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FIG.4-5 URVAKONDA TANK PARTICULARS 

(Stage V5 area and stage' -Vt Cumulative Volume) 
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CHAPTER 5 


Sin»!ARY OF RESOLTS, CCMICLUSIOHS AND SUGGESTIONS FOR FURTHER WORK 

6 . 1 General : 

The need for evaluation of groundwater balance is discussed 
earlier. Two methods for evaluation viz., the simple hydrologic 
budget method and the three dimensional finite difference method 
were presented in detail in the previous Chapter. In this Chapter 
the results obtained by the models are presented. The summary of 
results with conclusions and suggestions for further work are 
included in this Chapter. 

5.2 Results of the Hydrologic Budget Method of Analysis: 

The general equation for hydrologic budget method has already 
been explained in the previous Chapter. The two inputs to the 
basin are rainfall and seepage from canals. The outputs are 
consumptive use, evaporation from tanks and barren lands, outflows 
from the tanks. The balance between the input and the output 
gives the net change in storage. 

Using the available data for the years 1979 and 1980 for 
Havaligi basin, the above quantities are evaluated. Table 6.1 
shows the results. In this table, the rainfall volumes to the 
basin for the corresponding month is calculated using the average 
rainfall over the basin. The values of consumptive use are 
obtained by taking into account the various crops grown in the 
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basin. This includes evaporation from tanks, canals and from 
barren lands. The outflows from the tanks are available as input 
data. The q,uantity of seepage from canals is evaluated by using 
Darcy s Law in which a bed thickness of 30 cm, a coefficient of 
permeability of 0.4 m/month, and variable head depending upon the 
canal water level in that particular month are used. The 

algebraic sum of rainfall, consumptive use, outflows from tanks 
and seepage from canals is the net hydrologic budget. 

It is seen from the table that a total quantity of 50.75 
million cubic metres are added to the basin in the year 1979. 
Similarly, an amount equal to 32.68 million cubic metres are 
added. It may also be seen that in some months the net storage 
comes out to be negative. This means in those months this amount 
of water is removed from groundwater sotrage. 

5.3 Rainfall Recharge to the groundwater: 

One of the important component of groundwater budget analysis 
is rainfall recharge. The method of analysis to get the total 
recharge to the groundwater table has already been explained in 
Chapter 4. This total recharge includes rainfall recharge, canal 
recharge and recharge from other sources. In the present case, 
recharge only from two sources are considered namely rainfall 
recharge and recharge from canals. Table 4.8 gives the results of 
monthly rainfall recharge for the years 1979 and 1980, for each of 
the polygons into which the basin has been divided. The monthly 
values should be used rather cautiously as changes in soil 
moisture were not considered in computing these recharge values. 
In situations where soil moisture data are not available, it is 



63 


desirable to make an annual balance for groundwater recharge 
rather than monthly balance on the assumption that the soil 
moisture at the begining of each year will be the same. In our 
country, it is perhaps reasonable to assume that the soil moisture 
in the month of May (iust before the onset of monsoon) is the same 
in each year. By making such an annual balance for groundwater 
recharge, it is found that the rainfall recharge for the year 
1978-79 is 46.27 million cubic metres and for the year 1979-80 it 
is 40,82 million cubic metres respectively. Similarly the total 
hydrologic balance is also recalculated for the water years 
1978-79 and 1979-80 as 55.56 and 48.28 million cubic metres 
respectively. The difference between these numbers in each water 
year is because of seepage from canals which has already been 
discussed above. 

5.4 Results from Three Dimensional Analysis: 

The discretization of the basin has been explained in 
Chapter 4 . A total of 31 rows , 40 columns and 4 layers are used . 
Actually, last layer is a no-flow layer as required by the 
program. The aquifer characteristics for each of the cells is 
obtained by using the results of pump test conducted at a few 
selected points in the vicinity of the basin. 

The initial values of the heads used in the analysis are 
obtained from Fig. 4. 12. The boundary conditions for the basin are 
shown in Fig. 4. 11. 
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Using the above data, the three dimensional finite difference 
model was applied to the basin. The various types of results that 
we get from this model are: (i) the head distribution in different 
layers at different time periods (ii) the spatial distribution of 
drawdown at different time periods (iii) the overall volumetric 
Budget for each stress period. Since, the output is voluminous, it 
is not possible to give all the results. Table 5.2 gives a sample 
result of head distribution in layer 3 at the end of December 
1979. It is seen from the table that at some points the wells go 
dry; this is indicated in the table by a code number 1.0000E+30. 
The digits 9999. in the table represent inactive cells. Most of 
the tabular values are negative, as the datum was taken at the 
ground level. The cells in which the heads are more negative 
numerically are the ones which need attention in planning further 
pumping patterns. 

5.5 Conclusions: 

The objective of the present study was to implement the 
Modular Three Dimensional Finite Difference Groundwater Flow model 
on the mini computer namely Micro VAX II and on the DEC 1090 
systems. The second objective was to apply this model to a basin 
in Anantapur district of Andhra Pradesh. Implementation of the 
model has been done and the modifications that are required have 
been indicated. It is to be noted that, for successful 
application of the model to any basin, voluminous data are 
required. For example, the aquifer parameters are required not 
only in the lateral direction but also for different layers. Such 
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dei/silsd data may not be available for many basins. In this 
sense, there is a limitation of the applicability of the model. 
However, where such data are available, the model can be used for 
relaible results. It is further emphasised that the time varying 
data such as rainfall, runoff, evaporation, canal flows etc., are 
available for about ten years for reliable results. 

In the present case data for Havaligi basin was available for 
only two years. Hence the results of application of the model to 
the basin will only be approximate. 

In addition to the detailed three dimensional analysis, 
hydrologic budget method was also applied to the basin. This 
method will not take into account the variation in the hydrologic 
parameters of the basin and thus gives only the overall balance. 

5.6 Suggestions for Further Work: 

While applying the model, certain anamolies were noticed. For 
example, at some nodes, it was found that the head values were 
becoming positive, which indicates a spring. But, on the field no 
such springs were noticed. This may be because of paucity of the 
data. It is, therefore, suggested that in future work, in this 
basin, data over longer periods of time, say ten years, should be 
used. A part of the time series of the data may be used for the 
calibration and the remaining for validation. 
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